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Multiphonon processes in a model quantum dot (QD) containing two electronic states and sev- 
eral optical phonon modes are considered taking into account both intra- and inter-level terms. 
The Hamiltonian is exactly diagonalized including a finite number of multi-phonon processes, large 
enough as to guarantee that the result can be considered exact in the physically important region of 
energies. The physical properties are studied by calculating the electronic Green's function and the 
QD dielectric function. When both the intra- and inter-level interactions are included, the calcu- 
lated spectra allow for explanation of several previously published experimental results obtained for 
spherical and self-assembled QDs, such as enhanced 2LO phonon replica in absorption spectra and 
up-converted photoluminescence. An explicit calculation of the spectral line shape due to intra-level 
interaction with a continuum of acoustic phonons is presented, where the multi-phonon processes 
also are shown to be important. It is pointed out that such an interaction, under certain conditions, 
can lead to relaxation in the otherwise stationary polaron system. 
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I. INTRODUCTION 

The important role of phonons influencing the properties of systems based on semiconductor quantum dots (QDs) 
has been demonstrated in many works reviewed in Refsi^. The spatial confinement effect on optical phonons in 
QDs has been studied experimentally and understood theoretically 3 ' 45 6-7 . Acoustic phonons have also received 
some attention in connection with the low-frequency Raman scattering in spherical 8 ^ and self-assembled-^^ 1 QDs, 
and homogeneous broadening of the spectral lines 1 ^* 1 ^ 1 ^. Electron-phonon interaction in QDs remains, however, a 
controversial subject. While calculations performed for II- VI and III-V dots generally agree that the exciton-phonon 
(e-ph) coupling strength is reduced in nanocrystals compared to the bulk, they disagree regarding the numbers and 
trends in scaling of this interaction with QD size. (Electron-hole pairs will be traditionally called here excitons 
independently of the importance of the Coulomb interaction.) The calculated values vary substantially depending 
on the approximations used, but the Huang-Rhys parameter for the lowest exciton state usually does not exceed 
S 0.1 — 0.2 for II- VI for spherical QDs 1 ^* 1 ^ and is probably an order of magnitude smaller for III-V self-assembled 
dotsii. 

Turning to the experimental data, the exciton-optical-phonon coupling strength is most frequently obtained from 
PL spectra. Usually there is a large Stokes shift of the excitonic PL band with respect to the absorption in small II- VI 
QDs 118 , which can be explained by strong exciton-phonon coupling (Huang- Rhys parameter of the order of l^SLSi). 
However, as mentioned above, the calculated values happen to be one or two orders of magnitude smaller. Phonon 
replica in the PL spectra caused by the recombination of excitons in QDs were found to disagree with the well-known 
Franck-Condon progression both in their spectral position322*2I and relative intensitiesS ^*^ 3 ". Although both the 
strong Stokes shift and apparently large intensity of the first phonon satellite (relative to the zero phonon line) can 
be characteristic of the QD size distribution and not of a single doti&, these results indicate that the e-ph interaction 
in QDs must be considered carefully. Further evidence for this comes from the optical absorption measurements. 
LO-phonon related features were observed in absorption spectra of InAs/GaAs QDs and 2L02 4 *2£i2§i2 7 . and even 
3LC 28 phonon replica were found extraordinarily strong compared to the 1LO phonon satellite. At least some of 
these works were performed using single QD spectroscopy, so, ensemble effects can probably be ruled out in this case. 
Since no phonon replica were found in the corresponding emission spectra, there is no way to explain these results 
in terms of the Franck-Condon progression with any value of S. The temperature dependence of the homogeneous 
broadening of the absorption spectra of spherical II- VI QDs studied in Refs>i^S was found to have a clear contribution 
of optical phonons. This is unexpected, because coupling of optical phonons to a single electron or exciton level should 
lead to the appearance of discrete satellites and not to the level broadening. Thinking in terms of Fermi's Golden 
Rule, the broadening could be a lifetime effect owing to the electron (or exciton) transition to another state with 
emission or absorption of an optical phonon. However, it would require strict energy conservation in the electron- 
phonon scattering, i.e. exact resonance between the optical phonon energy and level spacing, which should be rather 
accidental. This kind of argument also justified the theoretical concept of "phonon bottleneck" , a very slow carrier 
relaxation which should be inherent to small QDs^. Nevertheless, an efficient phonon-mediated carrier relaxation 
has been reported in a number of works2i*22^a. All these experimental results imply that multiphonon processes are 
important and that the e-ph interaction in QDs must be treated in a non-perturbative way, even for the moderate 
values of the coupling constants coming out from the calculations^. An important ingredient to be included is the 
non-adiabaticity of this interactional 7 ^ leading to a phonon-mediated coupling of different electronic levels, even if 
they are separated by an energy quite different from the optical phonon energy. This is essential for understanding 
those experimental results which are in clear disagreement with the single level-generated Franck-Condon progression. 

Polaron effects in QDs have been studied theoretically in several recent papers22i2iiMi2£. The model considered 
in these works included two electronic levels and several Einstein phonon modes. The one-electron spectral function 
was obtained applying either self-consistent perturbation theory approximations 3 ^^ or exactly, using a combined 
analytical and numerical approach 3 ^. The results calculated using the perturbation theory approaches show shift 
and broadening of the levels, even for a sufficiently large detuning (defined as A' = A — Huq where A = (£2 — £1), 
£1 and £2 are the electron energy levels and huio is the phonon energy). However, as pointed out in Rcfs. 34-36 , this 
broadening is an artifact. The exact spectral function should consist of 5- functions (if no further effects are involved). 
Moreover, even if these 5-functions are broadened artificially, the self-consistent perturbation theory approximations 
are not able to reproduce the structure of the spectrum calculated exactly using the approach proposed in Rcf. 34 . 
The method of Refi^ based on the Gram-Schmidt orthogonalization procedure is, however, limited to the case when 
all the optical phonons have the same energy. In reality, confined optical phonon modes are characterized by different 
frequencies within the band of the corresponding bulk material. Usually there are a small number of such modes 
that interact more intensively with electrons or excitons^^ 7 .. The incorporation of this multiplicity seems to be 
important for comparison of calculated and experimental results. The limitation of the approach 3 " 5 * (based on the 
Davydov canonical transformation) consists in neglecting the multiphonon processes. 

Another important issue is the role of acoustic phonons. Even though some acoustic modes can be confined to the 
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QD and therefore have discrete energies, there should be a continuum of modes within a certain energy range. The 
interaction with this continuum should smooth out the polaron spectrum (otherwise consisting of S— functions, as 
noted above). The lineshape is not expected to be a simple Lorentzian for a QD, instead, acoustic phonon sidebands 
are formed for a discrete electronic leveli 4 ^. Considering the simultaneous interaction of confined excitons with 
acoustic and optical phonons can reveal some new effects and, to the best of our knowledge, has not been performed 
beyond the one-level modelSi. 

In this paper, we propose a non-perturbative approach to the calculation of the phonon effects on the electron 
spectral function, and emission and absorption spectra of a QD, based on direct numerical diagonalization of the 
Hamiltonian matrix including a small number of electronic levels and several optical phonon modes of different 
energy. Such a straightforward method was employed in Ref4& for phonons strongly coupled to an electron on a deep 
donor center. This approach is developed further here so as to allow for the incorporation of the interaction with 
(virtually all) acoustic phonon modes, under the condition that the inter- level coupling mediated by solely acoustic 
phonons can be neglected. We present calculated results which elucidate the influence of various relevant parameters, 
such as inter-level spacing, coupling strengths and temperature, on the optical spectra of a model QD. The paper 
is organized as follows. In Sec. II, we introduce the Hamiltonian matrix to be diagonalized in order to obtain the 
polaron states, and define the electron Green's function. In Sec. Ill, we derive the expression for the imaginary part 
of the exciton dielectric function and present the absorption and emission spectra calculated for different values of 
the relevant parameters. Sec. IV is devoted to the interaction with acoustic phonons and its incorporation into the 
calculation of the polaron spectra. We discuss the calculated results, make comparison with experimental data, and 
conclude in Sec. V. 



II. ELECTRON— OPTIC AL-PHONON INTERACTION: THE NON-PERTURBATIVE SOLUTION 

Our model system consists of two electronic levels, both coupled to N phonon modes (of frequencies u> u ). It is 
described by the Hamiltonian, 

2 N 2 N 

H l =Y J ^fa t + MtK + E £ 9tj a t a i (K + h) (1) 

i— 1 v— 1 ij u—1 

where af , a i are the fermion creation and annihilation operators for electrons (or holes) and , b v are the operators 
for phonons. The interaction with optical phonons occurs predominantly through the Frohlich-type mechanism and 
the Hamiltonian matrix elements between exciton states i and j are given by 4 ^, 

9ij = e J *|(r e , r h ) [(j> v {r h ) ~ <^(r e )] #j(f e , f h )dr e dr h (2) 

where § v is the electrostatic potential created by the z/-th phonon mode and ^ the exciton wavefunction of state i. 
In the following, we shall also use dimensionless interaction constants a"j — g^ / {hui v ) and omit the superscript when 
only one phonon mode is considered. 

If g\ 2 were equal to zero, the Hamiltonian would be exactly solvable, because there would be no interference 
between two electronic levels. The one-level model (known as independent boson model) has analytical solution 41 . 
The polaron spectrum, in the case of a single phonon mode, consists of equidistant peaks (separated by the phonon 
energy Tiujq). This is where the Franck-Condon progression comes from. 

The case of a\ 2 ^ 0, an = a 22 — f° r a single phonon mode, in the so called rotation wave (RW) approximation 4 ^ 
(which consists in neglecting the terms a^a^^ and a 2 a^b in (Q), also allows for analytical solution, with the 
spectrum given by 



£i l2 (m) = (A'/2) ± (A' 2 /a) + (m + l){g 12 ) 2 + (m + l)Huo (3) 

where m = 0, 1, 2, ... is the number of phonons in the mixed state. In addition to ©, there is a state with E = 0. 

The general case, which is of interest here, can be treated by mapping the many-body problem onto a single particle 
problem in a Fock space of higher dimensioniii 4 ^. It is natural to consider a basis \n\n2 {m^}) where n r = 0, 1 is the 
number of fermions on level r and m v is the number of phonons of mode v. In principle, the Hamiltonian matrix is 
infinite but, as it will be shown below, one can truncate the Fock space by allowing a certain maximum number of 
phonons for each mode and obtain a very accurate solution. Since we are interested in the case when there is a single 
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fermion in the dot, it is necessary to consider only the states |10 {m v }) and 1 01 {m„}). The necessary matrix elements 



are: 

JV 



SiS i: j + ^2 m v fiuj v 



N 

II 5 <ra v + 9ij Vm u + I5 m > umv+1 + Vr^ S mlm v -l] J[ 5 m>„ 



where i and j abbreviate the basis vectors 1 10 {m^}) (i — 1) and 1 01 {m u }) (i — 2). The dimension of the matrix is 
2toi • • • m v ■ ■ ■ rriN ■ For a small number of modes and a reasonable number of phonons for each mode, it can be easily 
diagonalized numerically. Given the eigenstates of the Hamiltonian matrix (denoted by |fc)), we can write down the 
electron Green's function^. In the canonical ensemble, 

GyW Z^f +G } E-(E k -E k ,)-ir, (4) 

where f3 — l/(fcsT) and Z = Tr exp(—/3H). Although rj in Eq.|@J should be infinitesimal, we suppose it to be a small 
quantity just for computational purposes. Equation (@J can be rewritten as 



r -0e. \ a t \ k ') ( k '\ gj | a 3 \ k ') ( fc, l 4 \ k ) 1 

UzA *> Z^ 6 \ E -(E k -E k ,)-ir, + E + (E k -E k ,)-ir,j 



(5) 



In the first term inside the brackets, the state has one electron and the intermediate state \k') has only phonons. 
The opposite occurs in the second term. Admitting that it costs an infinite energy to create a two-electron state, 
|11 {rriy}) cannot occur as the intermediate state in Eq.JSJ). Therefore, we can write for the diagonal elements of the 
Green's function: 

G U (E) =I^( e -^ + e-^,-^) \^H^})1 (6) 

kk ' E-(E k -J2 m v huj v ) - irj 

where ({m„}) are the eigenvectors expressed in terms of the basis vectors. From Eq.©, the fermion spectral 
density of states (SDS) is immediately obtained, 



P (E) = - Im £ G U (E) = 



i=l,2 i=l,2 

where pi is the partial SDS corresponding to the ?-th bare electronic level. 

The maximum number of phonons necessary to reproduce correctly several lowest-energy eigenvalues (which are 
of importance for relevant temperatures) depends on the coupling strengths but is not large. Taking just 10 as the 
maximum number of phonons allowed in the system, we obtained the eigenvalues coinciding with the analytical results 
of Ref4i (for an — , an = 0:22 = 0.2 — 0.3) and formula (for 012 = 0.2 — 0.3 , an = «22 = and using the 
RW approximation) to within 10 _3 meV. As an example, the SDS spectra calculated for the latter case (but beyond 
the RW approximation), are shown in Fig.l. 

III. EMISSION AND ABSORPTION SPECTRA 

One can calculate the optical absorption and emission spectra by using the Kubo formula for the frequency- 
dependent dielectric function^. The imaginary part of the dielectric function, which describes the absorption and 
emission properties, is related to the real part of the frequency-dependent conductivity by the relation 

47T 

Im e(uj) = —Re a(u>). (7) 
ui 

According to the Kubo formula, 

00 

Rea(c) = -^i- J dte^(j + (t)j(0)}. (8) 
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We are interested in describing the transitions from an electron state (labelled 0) in the valence band (which we 
assume is not coupled to phonons) to the states (i = 1,2) in the conduction band, and vice-versa. Thus, the current 
operator is defined as 

2 

3 = E Pot (af a o - a o a i) ■ ( 9 ) 
»=i 

The current-current correlation function is 

(j + (t)j(o)) = E PoiPoe <[4(*K(*) - 4(t)%(t)] - 4«o]>- (io) 

i,t'=l,2 

This expression can be developed into four expectation values. However, there are only two terms which contribute 
to the sum in Ea. (|10fl . This occurs because ajfc) vanishes if \k), an eigenstate of the Hamiltonian QJ, is a state with 
zero electrons. The non-zero terms are 

(4(t) ai (t)4a ) = ±J2^ Ek e lEkt (k lo+Oie-^o+Oo] k) 

k 

and 

{at{t)a Q {t)ata v ) = ± £ e _/JB *e <B ** (k \a+ a Q e~ im a+ a, \ k) 

k 

where a, = a,(0). These terms can be evaluated using the identity 1 = Yl \k')(k'\ by inserting the unity between the 

k> 

a operators. The result is: 

(i + (*)i(0)> = \ E E»e-^e J ^-^)* 

k,{m v } 



C$\{m v })C*{{m v }) - Cf ({m v })C*({m v }) 



(11) 



Using l|ll|) . performing the Fourier transformation in (jSJ) and substituting in Eq.(|7J), we get the following expression 
for the imaginary part of the dielectric function: 



Im e(u)) 



/ 2ire \ 
\mou) 



^ E {p*oiPoi' E e ^ m " %u3 »a 

' »,i' = l,2 fc,{m„} 



\{m v })CU{m v }) 



)(w - (E k - E ^)) +e~^ (fJ ^S, "^^^(o; + (S fe - £ m„fcj„))] } 



(12) 



where Z' = ^ exp — P^2 V rnvfuou . The terms in the second line of Ea. i|12f> correspond to the absorption and 
{m„} L J 

emission, respectively, of a photon of frequency u. Some absorption and emission spectra calculated for a hypothetical 
QD are presented in Figs. 2-5 (the parameters are indicated on the figures), that demonstrate the effects of the diagonal 
and off-diagonal coupling strength, inter-level spacing and temperature on the optical properties of the dot. In these 
examples, we considered the lower exciton level optically active and the upper one inactive (p2 = 0). Such a situation 
occurs in spherical II- VI (e.g. CdSe) QDs (lSelS'3/2 and ls e lP 3 / 2 states, respectively)^. 



IV. INTERACTION WITH ACOUSTIC PHONONS 



Interaction with confined longitudinal acoustic phonons, although mediated by a different (namely, deformation 
potential instead of Frohlich) mechanism, can be considered in the same way as that with optical phonons and should 
result in series of closely spaced by still isolated spectral peaks. However, for any kind of QDs embedded in a certain 
matrix, there must be a spectral region where acoustic phonons are essentially delocalized and their energy varies in 
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a continuous way. If the velocities of sound for the QD and matrix material are not very different, the density of 
acoustic phonon states should be similar to that of bulk crystals and these states can be characterized by a wavevector 
q. This is the case of self-assembled QDs^i and can be a reasonable approximation for a certain fraction of acoustic 
phonons in a spherical QD embedded in a dielectric matrix. In this section, we shall consider the interaction of an 
electron localized in the dot with acoustic phonons which are completely delocalized (see Appendix for details). The 
interaction, which occurs only inside the QD, is weak for each phonon mode (since it contains a factor Vqd/V, V is 
the volume of the whole system). However, since the number of modes is virtually infinite, perturbation theory may 
fail. We shall avoid it. For this, we will have to neglect coupling between different electronic levels through acoustic 
phonons and consider a single electronic level coupled to arbitrary number of phonon modes. This approximation 
corresponds to the independent boson model^i, which is normally considered for optical phonons (see Sec. III). 
Recentlyi^S, it was applied to acoustic phonons in a QD. Based on the exact solution of this model, we shall propose 
a different procedure which will lead us directly to the self-energy function of the electron or exciton. Later, this 
self-energy will be re-interpreted for the optical-phonon polaron. 

The Hamiltonian of a system consisting of one electronic level and a continuum of acoustic phonon modes is: 



H 2 = a + a(e + ^2g<i(b+ + b^j + ^huj^b+b^ 



q q 

41 



It can be diagonalized by the transformation to new bosonic operator 

£? q = 6 q + a q a + a 
where a q = (<? q /?kj q ). The energy spectrum is given by 

S({iV q })=4 + ^na; q iV q 



(13) 



(14) 



(15) 



with e' q = s — J^ q a^liWq. Eigenstates of the Hamiltonian 113|l can be expressed in terms of the pure phonon states, 

W = E^K). ( 16 ) 



Acting with the operator B+B on the wavefunction (|16H we can find recurrent relations for the coefficients C 



N ■ 



N = (m + al)(C%) 2 + a c 



N 



In linear (in the interaction with a single phonon mode q) approximation, we get from l|16|) and l|17fl : 



W = PN, V^q«ql^q ~ + \ N l) + V N 1 + ^ql^q + 1) 



where 



(17) 



(18) 



PN„ = 



'1 + (2A r q + l)a\ 

The one-electron Green's function corresponding to the Hamiltonian (|13|) can now be calculated using Eq. , 



G{E) 



E 

{m q } 



]l{<WiV q + «q[^q<Ww q -l + (Nq + 1 )<5m q , JV q + l] } 



E - e' Q - Y, ftUq(Nq - TO q ) - ITj 

q 



xexp [-^(2A q + l)^ 



(19) 



where the angular brackets stand for thermodynamical average and the exponential factor appeared from the product 
]^[pAr q . Taking into account only one-phonon processes, the thermodynamical average approximately replaces iV q -s 
q 

with the corresponding Bose factors A q = ^exp (/3hujq) — lj and the Green's function can be written as: 



G{E) = { 



n, 



«q^q 



« q (iV q + 1) 



E — £q — 17] Z -— ' IE — E'q — Hujq — IT) E — Eq + Huiq — IT\ 



} 
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+ (2^ + iKj . (20) 
q 

However, this equation (|20|l can be used only for quite a weak interaction or at very low temperatures. Although 
the interaction is weak for each phonon mode (all a q are small), the number of modes is large and the effective 
electron-phonon interaction is strong. A typical number of phonons interacting with the electron can be estimated as 
Q = (2Aq + l)a!q and, for the model presented in Appendix, Q » 1 for temperatures usually used in the experiments. 
In these conditions, it is necessary to use the general equation (|19|1 . 

The evaluation of the Green's function from Ea. H19l) can be made using a Monte-Carlo procedure. Instead of 

summing over all 3 configurations {fTi q } (N is the number of acoustic phonon modes), one can generate some N^ic 
most probable configurations with m q = (A q — 1), Aq, (A q + 1), occurring with the probabilities a q A q p^ , p 2 ^ , 

and a q (A q + l)p% , respectively. Averaging over a sufficiently large Nmc of such configurations, G(E) is obtained 
including all the many-phonon processes. 

Let us define the electron self-energy taking into account the interaction with acoustic phonons, 

Y,(E)=E-G- 1 {E). (21) 

In the one-phonon approximation, the explicit expression for the self-energy is: 

E(£) =4 -(£-4)5> q (2AT q + l) 



t r a q A q a q (Aq + 1) 1 

+ (E-e ) ^ y- _ - - + - - - - j . (22) 



Tq — hojq — vt) hj — e + hojq — irj 

The spectral dependence of the self-energy corresponding to the Green's function (Eq.J^OJl) obtained using the Monte 
Carlo procedure is presented in Figs. 6, 7. For comparison, we also calculated the self-energy using Ea. (|2*2*|l . The 
interaction constants a q derived in the Appendix were used in these calculations and the material parameters were 
taken as for CdSe, except for the deformation potential constant a c , which was taken about 2 times smaller than the 
bulk value of the relative volume deformation potential between the valence and conduction bands^t Such a choice 
is justified by the fact that only a fraction of acoustic phonons can be described by a propagating wave assumed in 
the Appendix. 

Assuming that the broadening of the electronic levels produced by the interaction with acoustic phonons is small 
compared to the electronic level spacing, it is possible to include this interaction in our scheme of consideration of 
polaronic states proposed in the previous sections. Let us consider now the full Hamiltonian, which includes the 
interactions with both optical and acoustic phonons, 

2 2 N 



H * = E + E 4° K + & q) a t a i + E E f« tfs 



/V 



+ Y t ^b+b^ + Y l M+K (23) 

q v— 1 

where a q ^ denotes the acoustic phonon coupling constant for electronic level i. The Hamiltonian (|23|l can be rewritten 
in terms of the polaron states \k) by introducing the corresponding (fermionic) annihilation and creation operators 
A^,A k . The electron-LO-phonon Hamiltonian is then reduced to 

H x = Y,E k A+A k 

k 

Using the expansion of the bare electronic states in terms of the polaron ones, 

i^EE^^}))*^' 
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the term in (12,'j|) representing the interaction with acoustic phonons can be written as 

E E^^^E E (pi ({m„}M+A„ 

»=1,2 q fc.fc' {m„} 

= E 5 q (^ + 5 q) E A i A k + E A 5 q (^ + fe q ) E ^.fe'^fc^fc' ( 24 ) 

q k q fe,fe' 

where g q = (pq 15 + ,g q 2) )/2, A 5q = ( 5q 2) - g q x) ), and 

tkM = E {(^ (W))* C i' (M) ~ ( C * (W))* (M) }■ 

{m„} 

Neglecting the last term in l|24f) . that is, assuming that the coupling to acoustic phonons is approximately the same 
for both electronic levels, we come to the Hamiltonian 

H 3 = E A t A k (p h + E 5q(^ + bj) + E t*>«bp v (25) 
k q q 

It can be diagonalized exactly in the same way as . The poles of the electron Green's function © (as well as the 
absorption and emission peaks) are shifted according to the replacement Ek — > E' k , 

E' k = S fe - E "q^q + (E W-Nq - TO q))> 
q q 

where, as before, (7V q — m q ) can be considered a random variable taking the values 1, and -1 with the probabilities 
a^Nqpjy , p 2 jj , and a q (iV q + l)p% , respectively. The result is a similar broadening and a downward shift of all the 
polaron states contributing to the spectral density of states, absorption and emission, just like it happens in the one- 
level independent boson models. It is equivalent to ascribing a spectral variable-dependent self-energy, Y,(E — E p ) 
(given by Eas. i|21f> or (|22|l % ). to each pole of ©. As before, we can re-interpret this result for exciton-polaron, under 
the condition that we have no more than one exciton per QD. An example of spectra showing this effect is presented 
in Fig. 8. The non-diagonal term with tk,k' leads to acoustic phonon-mediated mixing of the polaron states. It will 
be considered in a future work. 



V. DISCUSSION 



Let us start the discussion by emphasizing that only in a hypothetical case when the phonon-mediated coupling 
of the lowest energy exciton state to all other states can be neglected (for example, in the limit of extremely strong 
confinement, such that A >> Hujq, or if ot\i is small because of the symmetry of the corresponding wavefunctions), 
one can expect to observe Franck-Condon progressions in the emission and absorption spectra associated with this 
stated. This is the only case when the Huang-Rhys parameter describes the spectra in entirety. As it is obvious 
from Eq. @ , the diagonal coupling is proportional to the (integrated) difference between the electron and hole charge 
densities, which is non-zero mostly due to the fact that the conduction and valence bands in II- VI and III-V materials 
have different symmetry. Even if some further effects (like the presence of defects or strain, for example) contribute 
to the separation of the electron and hole clouds in spaced, one can hardly expect an order of magnitude increase of 
this parameter, compared to the calculated values cited in the Introduction. It means that the intensity of the LO 
phonon sidebands must be rather small and monotonically decreasing with the increase of the number of phonons. 

However, in most cases the phonon-mediated interlevel coupling, at least between the subsequent exciton states, 
should be important. The off-diagonal interaction constant is likely to exceed the diagonal ones. For instance, if the 
two exciton levels arc different by the hole state (e.g., ls e lS 3 / 2 and ls e lP 3 / 2 states in a spherical II- VI QD), it is 
easy to see from Eq.@ that (neglecting the Coulomb interaction between the electron and hole), 

9X2 = ej ^* lh (r h )Mrh)^2h(r h )dr h , (26) 

that is, there is no compensation effect for this interaction. The spectra calculated assuming only the off-diagonal 
interaction (Figs. 1,2) show characteristic features known as Rabi splitting, which is obtained already in the RW 
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approximation (see Eq.J3|)) and has been observed experimentally 31 . Contrary to the RW approximation, the exact 
results also show a downwards shift of the spectral lines (see Fig.l. In the RW approximation, the most intense peak 
should be situated exactly at E — Eo = 0). Even if the upper level is dipole forbidden, its presence manifests itself by 
the spectral features seen in Fig. 2. Note that the off-diagonal coupling is not resonant as one might expect thinking 
in terms of the Fermi's Golden Rule. The manifold produced by the coupling persists even for quite a large detuning 
A' > TujJq. Curiously, it is reminiscent of sidebands which might appear due to interaction with confined acoustic 
phonons. 

When both diagonal and off-diagonal interactions are present, the effect is not just a sum of those owing to each 
of them, but also some additional spectral features appear (see Figs. 3-5 showing the effect of different parameters). 
The stronger are the interactions, the richer the structure of the spectra. It should also be noted that the symmetry 
between the absorption and emission spectra, which is characteristic of the one-level independent boson model (see the 
upper panel of Fig. 5), disappears when the off-diagonal interaction is included. This holds also at low temperatures 
where a small number of spectral features are seen (Fig. 3, lower panel). Those below Eq in the emission spectrum are 
due to the diagonal coupling to the lowest energy exciton state. The one seen in the absorption originates from the 
inter-level coupling. 

Let us turn to the acoustic phonons. The explicit calculation leads to a self-energy which depends on the spectral 
variable in a sophisticated way. Accordingly, the line shape is quite different from the simple LorentzianMi^. The 
noise seen in the calculated spectra of the self-energy (Figs. 6, 7) is due to the Monte Carlo procedure used when 
each configuration {?n q } of phonons produces a (5-function and the number of Monte Carlo runs is finite (Nmc — 
20000 for Figs. 6, 7). The imaginary part of the self-energy vanishes at e' Q , so, the zero-phonon line (ZPL) at low 
temperatures remains a sharp 5-like resonance situated between two asymmetric phonon sidebands, in accordance 
with experimentpiSii 3 . and previous theoretical studies^*^. Only for high temperatures the SDS generated by a single 
electronic level takes a bell-like shape similar to the Lorentzian. The calculated results also confirm the statement 
concerning the importance of the multiphonon processes which increase in smaller QDs (Fig. 6). Turning to the 
exciton-polaron spectrum (Fig. 8), the acoustic phonon-related broadening smears the discrete structure associated 
with optical phonons The multiplicity of the confined optical phonon modes also helps this. Below, we shall explain 
some previously published experimental findings, which were not clearly understood before, in terms of our calculated 
results. 

"Strange" phonon replica 

An experimental evidence of spectral features whose distance from ZPL is smaller than the LO phonon energy was 
found in several works22i2ii2i2&. Thinking in terms of the Franck-Condon progression, such peaks were attributed to 
additional (interface or disorder-activated acoustic) phonon modes, for some reason strongly coupled to the exciton^. 
As it can be seen from Figs. 3-5, in the presence of inter-level coupling, there are many spectral features that are not 
separated from ZPL (the most intense peak in all spectra) by (a multiple of) the optical phonon energy. Some of them 
persist at low temperatures. Such "strange" replica can be generated by a single phonon mode (long-wavelength LO 
phonon, neglecting its confinement) and their spectral positions are determined by the e-ph coupling constants. Thus, 
there is no need to invoke extra modes (for which it would be hard to justify their strong coupling to the exciton) 
in order to explain the discussed spectral features. The phonon-mediated inter-level coupling provides a simpler and 
more plausible explanation. 

Anomalously strong 2LO phonon satellite 

As mentioned in the Introduction, several groups observed apparently phonon-related features in the absorption 
(but not in the emission) spectra of resonantly excited self-assembled QDs^2£*2£i2L2£. Given the relative weakness of 
the exciton-phonon interaction in such dots (S << 1), the most difficult to explain was the fact that the most intense 
satellite was not the 1LO phonon sideband but rather the 2LO or 3LO one. This effect was qualitatively understood 
as being due to a resonant coupling of the corresponding phonon replica of the ground state to the excited state of 
the excitonS&. Considering a QD ensemble, one can think in terms of a multiphonon "filtering" of inhomogeneously 
distributed excited states^. However, when a single QD spectroscopy is usedSi, such a resonance of the inter-level 
spacing with a multiple of the phonon energy is rather improbable. (Even though there are several confined optical 
phonon modes with slightly different energies, capable of considerable coupling to the exciton, this dispersion is only 
of the order of lmeV). The assumption of a reasonably strong off-diagonal coupling of the lowest and higher exciton 
states (not in a close resonance with a certain number of phonons) can account for anomalously strong nLO phonon 
sidebands. Fig. 9 shows how peaks separated by energies approximately equal to that of the optical phonon can appear 
in the absorption spectra. We took a value of A=55meV, typical of InAs/GaAs self-assembled QDs, and considered 
two optical phonon modes with the energies of 32 and 30 meV in order to simulate the experimental result of Lemaitre 
et aS&. (The two modes can be interpreted as LO phonon in pure InAs and InAs-like LO phonon in the InGaAs 
alloy near the QD boundary, respectively). The values of a" 2 — 0.2 and 0.15 used in this calculation do not look 
extraordinarily high taking into account the argument presented above (see Eq. l|2t)|) ) and the possible contribution 
of the optical deformation potential mechanism in this interaction. As it can be seen from Fig. 9, the second strong 
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absorption peak is not separated from ZPL by the energy A and can be called "2LO satellite" , since the weaker 1LO 
and 3LO ones also show up in the spectrum. The calculated spectrum is in a good agreement with the experimental 
observation of Ref— . 

Homogeneous broadening of absorption lines 

In structures of higher dimensionality, the temperature dependence of the homogeneous linewidth is described by 
the following relation, 

r(T) - T(0) = lac n ac (T) + lL on L o{T) (27) 

where 7 ac and ^lo are some constants and n ac and n^o the Bose factors for the acoustic and optical phonons, 
respectively. There is no reason for Ea. (|27f) be a good approximation for QDs, because it relies on the existence of a 
continuum of electronic states and also is limited to one-phonon scattering processes. Still it was used i n 15 i 29 and some 
other publications to qualitatively describe the dependence T(T)) extracted from experimental absorption spectra of 
spherical QDs. It was found that this dependence, in the region from 77 to 300K, is much stronger than linear, similar 
to that predicted by the second term in Ea. l(77|) . suggesting that there is some broadening related to the optical 
phonons. As it has been pointed out in the Introduction, interaction of localized electrons with optical phonons 
produces only additional discrete spectral lines. Consequently, the optical phonon related part of the broadening 
should not be understood literally, strictly speaking, the broadening occurs only because of the continuum of acoustic 
phonons. Our calculated Fig. 8, where the interaction constants were taken as for a typical spherical CdSe QD studied 
in Ref»2^, demonstrates this effect. The interaction with several optical phonon modes results in a series of closely 
spaced discrete peaks (see inset in Fig. 8)) which is camouflaged by the acoustic phonons. The effective homogeneous 
broadening of the spectral lines may apparently depend on the equilibrium number of acoustic and optical phonons. 
Nevertheless, Ea. (|27[l is too simple to describe this effect even qualitatively, as found in Ref— using single QD 
spectroscopy. 

Owing to the apparent similarity with systems of higher dimensionality, some authors associate the broadening 
discussed above with the optical phonon-assisted carrier relaxation. However, as it has been pointed oui£2£&, this 
is not a lifetime (or irreversible scattering) effect. Indeed, the polaron states considered here are stationary states. 
The phonon-assisted relaxation can probably occur through an anharmonic decay of the phonons participating in the 
formation of the polaronic states, as suggested in Refs— and calculated in Ref— . Our consideration of the polaron 
spectrum in the presence of acoustic phonons shows, however, another possibility, namely, the acoustic phonon- 
mediated transitions between different polaron states owing to the non-diagonal term in (|24H . The next step should 
consist of analyzing the efficiency of this mechanism. It is worth noting at this point that, contrary to the opinion of 
the authors of Ref— , such acoustic phonon-mediated transitions in the polaron spectrum should not be subject to 
phonon bottleneck, because there are plenty of polaron states due to different optical phonon modes differently coupled 
to the exciton. Many of them (which not necessarily show up in the spectra corresponding to thermal equilibrium) 
are separated by small energies within the band of acoustic phonons efficiently interacting with the polaron. 

Up-converted PL 

The up-converted or anti-Stokes photolumincscence (ASPL) is the emission of photons with energies higher than 
the excitation energy (E exc ). This effect was observed for colloidal II- VI QDs and discussed in several recent 
publications^iS^i. Additional information can be found in Ref. 55 . The ASPL occurs when an ensemble of QDs 
is excited at the very edge of the absorption spectrum, below the normal (i.e. excited with high energy photons) PL 
band. The principal experimental facts concerning this effect are the following: 

(i) The ASPL intensity increases linearly with the excitation power—, which can be rather low. 

(ii) The blue (anti-Stokes) shift between the ASPL peak and E exc does not depend significantly upon the QD size 
if E exc is chosen proportionally to the absorption peak energy (which depends upon the size)— However, the shift 
increases with temperature and can range approximately from 20 to 150me\— 

(iii) If E exc increases (approaching the absorption peak) the ASPL also moves continuously towards higher 
energies 53 . Its intensity increases, and finally the spectrum transforms into the normal PL band—. 

(iv) The ASPL intensity increases strongly with temperature. 

Possible ASPL mechanisms were discussed in Refs— but no definite conclusion was made. According to (i), 
processes like two-photon absorption and (since the possibility of emergence of more than one exciton per QD is 
negligible) Auger excitation can be excluded in this case. Therefore, it was suggested that incident photons excite 
electrons to some intermediate sub-gap states from which they eventually proceed to the higher energy (luminescent) 
states through the thermal effect. Some obscure surface states (SS's) for which, as admitted by the authors of Ref— , 
there is no any direct evidence, were suggested as those responsible for the sub-gap absorptio n 53 ! 54 . 

From our point of view, the participation of sub-gap SS's (even if they exist in the passivated colloidal QDs) is 
unlikely in virtue of the experimental results (ii) and (iii). In fact, it would require the SS's energy have approximately 
the same dependence on the QD size as that of the confined electronic states in the dots, in obvious contradiction with 
what should be expected from the general point of view. At the same time, there exist naturally formed sub-gap states, 
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which are separated from the fundamental absorption line (i.e., ZPL) by energies which are only weakly dependent on 
the QD size. These are the red-shifted optical phonon replica. Fig. 10a presents the lower-energy side of a calculated 
QD absorption spectrum showing two sub-gap bands (designated as "-1LO" and "-2LO") through which the dot can 
be excited. The excited QD then will emit a photon most likely having the ZPL energy. The probability of such 
an up-conversion process increases with temperature because so does the integrated intensity of "-1LO" and "-2LO" 
absorption bands (shown calculated in Fig. 10b). The experimental temperature dependence of the ASPL intensity 
measured in RefS can be understood taking into account that, at a certain temperature, further red-shifted satellites 
(whose intensity depends more strongly on the temperature) become more efficient. The situation is complicated by 
the distribution of the QD size, so that E exc can match different "— nLO" bands of dots of several different sizes. This 
agrees qualitatively with the experimentally observed increase in the anti-Stokes shift with temperature^. Modeling 
of the ensemble effects involved in the ASPL phenomenon is beyond the scope of this paper and will be considered 
in a future work. Nevertheless, we believe that our model explains, at least qualitatively, the principal experimental 
facts concerning the ASPL. 

In conclusion, we proposed a non-perturbative approach for the calculation of the polaronic effects in QDs, which 
allows the consideration of electronic levels coupled through the interaction with several confined phonon modes. Using 
this approach, we were able to show that the polaronic effects are significant even when the interlevel spacing is quite 
away from resonance with the optical phonon energies. We demonstrated that this opens the possibility to account for 
some previously published and not clearly understood experimental data. We also presented an explicit calculation 
of the spectral line shapes because of the diagonal interaction of acoustic phonons with the original electronic levels 
and suggest that the same interaction may be responsible for relaxation of the polaron entity to lower energies. 
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APPENDIX: ELECTRON- ACOUSTIC-PHONON COUPLING RATES 



The following simple model was used to estimate the rates a q of electron coupling to longitudinal acoustic phonons. 
Assuming that LA phonons arc completely dclocalized, that is, neglecting acoustic impedance at the interface between 
the QD and matrix, the displacement can be written as 

Uz = ( V /2 Re e *(g*-<" fl t) 

z \2uj q poV J 

where po is the density. For the deformation potential interaction mechanism, the dimensionless coupling constant 



« q = £^ J |Vfdivu(q)dU, 



where a c is the bulk deformation potential and ip the electron wavefunction. Considering the lowest conduction band 
state in a spherical QD of radius R, with infinite barriers, the wavefunction is^fe 



1 sin(7rr/i?) , , 



Using these expressions, we obtain 

1 1 

. 2~huj q poVc 



"q-"c( ot ,, „ jjrs ) I(qR) (A.l) 



where q is the longitudinal sound velocity, 



I(y) = 2tt 2 f j (y ■ x)jl{nx)x 2 dx (A.2) 
Jo 



and jo is the spherical Bessel function. The integral in (|A.2fl can be expressed in terms of the integral sinus as 

i{Si( 2/ )-i[Si(y + 2 7 r) + Si(|y-2 7 r|)]}. 
The function I(y) decreases from 1 for y = to nearly zero for y > 5. 
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FIG. 1: Calculated partial spectral densities of states corresponding to two electronic levels separated by A = 30meV 
interacting with one optical phonon mode with huo = 30meV (that is, A' = 0). Only non-diagonal coupling was included 
(ai2 = 0.3). Note the Rabi splitting in the polaron states and their shift to the lower energies. 
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FIG. 2: Calculated absorption spectra for different values of level spacing as indicated. All the other parameters are the same 
as for Fig.l. 
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FIG. 3: Absorption (left) and emission (right) spectra for a system of two exciton levels separated by A = 30meV calculated for 
different temperatures. The excitons interact with two optical phonon modes of the energies Tiuji — 25meV and %u)2 = 25.2meV . 
The electron-phonon coupling constants are indicated on the figure. The origin of the energy axis is chosen at the energy of 
the lowest exciton state (Eo). 
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FIG. 4: Absorption (left) and emission (right) spectra calculated for different (compared to Fig.3) values of the diagonal 
coupling constants. 
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FIG. 5: Absorption (left) and emission (right) spectra calculated for different values of the non-diagonal coupling constant 
(see the explanation to Fig.3). 
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FIG. 6: Imaginary (left) and real (right) parts of the self-energy for a single electronic level interacting with acoustic phonons, 
calculated for different temperatures using the Monte-Carlo technique explained in the text. The dashed curve was plotted 
using the one-electron approximation (Eq.(22)). The coupling constants were calculated as explained in the Appendix. The 
necessary material parameters were taken as follows: po = 5.8g/cm 3 , ci = 3.8km/s and a c = 3.5eV. The QD radius i?=2nm. 
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FIG. 7: Imaginary part of the self-energy calculated for a i?=4nm QD as explained for Fig. 6. The dashed curve was plotted 
using one-electron approximation (Eq.(22)). 
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FIG. 8: Absorption spectra calculated for an _R=2nm CdSe QD including both optical and acoustic phonons for three different 
temperatures. Two optical phonon modes of the energies Tiuji = 25meV and Tiui2 = 25.8meV and dimensionless interaction 
constants (divided by the corresponding phonon energies) a\ 1 = —0.015, Q22 = 0.027, a\ 2 = 0.295; afx = 0.063, ol\ 2 — 0.053, 
«i2 = 0.095 calculated in Reft— were used. The acoustic phonon parameters were taken as for Fig. 6. The level spacing of 
A = 80meV was assumed. The inset shows the T — 300K spectrum without acoustic phonon contribution (at T = 10K only 
the main peak is seen). 
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FIG. 9: Low temperature absorption spectrum calculated for an InAs QD with the level spacing of A = 55meV assuming 
two optical phonon modes (the energies are given on the figure) and neglecting coupling to acoustic phonons. The interaction 
constants used in this calculation are: a\x = 0.005, ot\ 2 = 0.005, a\ 2 = 0.2; afj = 0.005, Q22 — 0.005, »i 2 = 0.15. Notice that 
three features above the zero-phonon line are separated by nearly the same energy, approximately equal to 31meV. The inset 
shows the corresponding emission spectrum. 
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FIG. 10: (a) Low-energy part of the absorption and emission spectra calculated for a hypothetical CdTe QD considering three 
optical phonon modes with parameters given on the figure. The acoustic phonon parameters were taken as for Fig. 6. The level 
spacing A = WOmeV , temperature 300K. (b) Temperature dependence of the integrated intensity of two sub-gap bands in the 
absorption spectrum (a) (lines) and experimental data of Reft— ■ (points) showing the temperature dependence of the ASPL 
peak amplitude. 



